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Abstract. The processes described in the title always have reversible stationary distributions. 
In this paper, we give sufficient conditions for the existence of, and for the nonexistence of, non- 
' reversible stationary distributions. In the case of an i.i.d. environment, these combine to give 

a necessary and sufficient condition for the existence of nonreversible stationary distributions. 



j> ■ 1. Introduction 

o : 

Transport phenomena for noninteracting particles in one dimensional environments is a 
well studied subject in the contexts of classical and quantum systems. However, the influence 
of interactions among particles in these situations is considerably less well understood. In 
this paper, we consider one of the more important models of particle motion with interaction 
- the exclusion process. In particular, we will study the exclusion process on Z 1 with 
nearest neighbor jumps with probabilities pi G (0, 1) and qi = 1 — Pi from i to i + 1 and 
% — 1 respectively. This is the continuous time Markov process r\ t on {0, 1} Z with formal 
generator 



O 

c3 



X 



£f(v) = E {»?(*) [1 - V(i + l)]pt + V(i + 1)[1 - »7(0]ft+i}[/fai,i+i) - f(v)i 

i 

where ?7i,i+i is obtained from 77 by interchanging the ith and (i + l)st coordinates. The 
exclusion process has been the object of a lot of attention over the past 35 years, primarily 
in case the transition probabilities are translation invariant. Here, as the title indicates, we 
investigate problems where the pi are inhomogeneous, with particular applications to the 
case in which the pi are i.i.d. random variables. Among the few rigorous papers dealing 
with spatially inhomogeneous asymmetric exclusion processes are Liggett (1976) and Jung 
(2003). 
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The i.i.d. model was investigated in a certain approximate fashion in Tripathy and Barma 
(1998) and, more recently, in Harris and Stinchcomb (2004). In general, these authors have 
studied the "phase diagram" of (maximal) current flow as a function of an equilibrium 
particle density parameter in the presence of disorder. Certain interesting phenomena have 
been uncovered in these works. In particular, a symmetric flat region in the above mentioned 
response curves, indicating a forbidden interval of densities, and the observation, primarily 
numerical, in the earlier reference that in the presence of a random locally preferred direction 
of flow, the current density vanishes with increasing system size. Our results bolster some 
of these conclusions: We demonstrate the existence of current carrying states whenever, 
for some e > 0, po > | + e. Moreover, the flux in these states vanishes linearly with 
e. Furthermore, we vindicate conclusively the later phenomenon. In particular, whenever 
Po < qo and po > qo have positive probability - or even if po = \ is in the support of the 
disorder distribution - we show that the current indeed vanishes as the system size tends to 
infinity. Furthermore, under the general condition of zero current, we can characterise all 
the invariant measures. 

This picture is in sharp contrast to the non-interacting version of this problem: At 
the end of this paper, we will consider briefly the non-interacting case, and show that 
there is a stationary distribution for the system with nonzero flux whenever the pj's are 
i.i.d. and Elog(pi/qi) exists and is nonzero. This difference between the interacting and 
noninteracting systems might be a bit surprising, since at least at low densities, one might 
expect the two systems to have similar properties. 

Returning to the exclusion process, we first note that reversible stationary distributions 
always exist in our situation. To define them, let 7Tj be defined by taking -kq > and then 
T^iVi = ^i+iQi+i f° r au i £ Z 1 • The corresponding reversible measure is the product measure 
u a , where 



(See Section VIII. 2 of Liggett (1985), for example.) According to Theorem 2.1 of Jung 
(2003) , these are extremal in the class X of all stationary distributions if and only if 



When this sum is finite, extremal stationary distributions /z n (with n > or — oo < n < oo, 
depending on the situation) are constructed from these product measures in the following 
way: //„(•) = v a (- | A n ), where 



(1.1) 




(1.2) 





where in the third case, T is chosen so that ^2, ieT o.i < oo and Yligri^ ~ ai ) < 00 • These 
conditional measures do not depend on a, since an irreducible positive recurrent Markov 
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chain has a unique stationary distribution. In the third case, {p n } does not depend on the 
choice of T, except for a possible relabelling. 

In the spatially homogeneous case pi = p, the extremal stationary distributions are 
completely known (Liggett (1976)): 

X e = {z/ p , < p < 1} 

if p = 1/2 and 

X e = {v p , < p < 1} U {p. n , — oo < n < oo} 

if p 7^ 1/2. Note that in the latter case, v p is not reversible. (Here v p denotes the homoge- 
neous product measure of density p.) 

The spatially inhomogeneous case in which Z 1 is replaced by {0, 1, 2, ...} was treated in 
Liggett (1976). In that case, the result is that all stationary distributions are reversible. 

Our main objective in this paper is to say what we can about the following question: 
In the spatially inhomogeneous case on Z 1 , when do nonreversible stationary distributions 
exist? We will give sufficient conditions for the existence and for the nonexistence of such 
distributions; they become necessary and sufficient in the case of i.i.d. p^s. 

An important tool in discussing this issue (as in many involving the exclusion process) is 
the flux. For p G X, this is defined by 

(1.3) <p(p) = Pip{r] : r](i) = 1, r](i + 1) = 0} - q i+ ip{v ■ = 0, n(i + 1) = 1}. 

The fact that this quantity is independent of i can be checked by using J Cfdp = for 
f(v) = vii)- m case Pi = Pi the flux under v p is (p — q)p(l — p), for example. 

It is easy to check that the flux is zero for the reversible stationary distributions described 
above. Our first result is a converse to this observation. 

Theorem 1. Suppose p G X e and 4>{p) =0. If p is not the pointmass on n = or on 
T] = 1, then p = v a for some a satisfying (1.1) if (1-2) is satisfied and p = p n for some n 
if not. 

For applications of Theorem 1, it is useful to make the following observation: 

Proposition 1. Take e > 0. Suppose that for each n there exists a k so that pi > | — e for 
all k < i < k + n. Then 4>(p) > — e for all p G X. Similarly, if for each n there exists a k 
so that pi < \ + e for all k < i < k + n, then 4>{p) < e for all p G X. 

Theorem 1 and Proposition 1 will be proved in Section 2. Combining these two results, 
we obtain a sufficient condition for all stationary distributions to be reversible: 

Corollary 1. Suppose that for every e > and every positive integer n there exist k and I 
so that pi > \ — e for all k < i < k + n and pi < | + e for all I < % < I + n. Then p G X 
implies that p is reversible. 

This conclusion that all stationary distributions are reversible under the much stronger 
assumption that 

lim pi = — 

i— >±oo 2 
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is a consequence of Theorem 1.1 of Jung (2003). This conclusion also follows from Theorem 
1.2 in that paper (assuming inf iPi > 0, infj qi > 0) in all cases other than (a) lim^-oo %i = 
0, limi^ +00 TTi = oo or (b) lim^-oo %i = oo, limj_ >+00 = 0. When the p^s are i.i.d., these 
excluded cases are of course the prevalent ones: case (a) occurs when E\og(po/qo) > and 
case (b) when E\og(qo/po) > 0. 

To prepare for the next result, we define the exclusion process on [m,n] = {m,...,n} 
with boundary conditions by allowing the usual transitions in {m, ...,n} together with the 
transitions — > 1 at site m at rate p m -i and 1 — > at site n at rate p n . This is a finite 
state irreducible Markov chain, and hence has a unique stationary distribution /U m n . 

Theorem 2. TTie 0(// m ,n) 

is o?7. increasing function of p m —i, p m , ...,p n . 

Combining this with the known behavior of \im,n i n the homogeneous case yields the 
following sufficient condition for the existence of a nonreversible stationary distribution: 

Corollary 2. If for some e > 0, pi > \ + e for all i, then there exists a (nonreversible) 
|i£l with 4>(n) > |. 

This answers a question raised near the end of the introduction to Jung (2003). Theorem 
2 and Corollary 2 will be proved in Section 3. 

Combining Corollaries 1 and 2, we have the following result for the exclusion process in 
which the p^s are chosen randomly in an i.i.d. fashion: 

Theorem 3. Consider an exclusion process with i.i.d. pi 's. The following hold with prob- 
ability 1: 

(a) If for every e > 0, P(po > \ — e) > and P(po < \ + e) > 0, then all stationary 
distributions are reversible. 

(b) If for some e > 0, P(po > \ + e) = 1 or P(po < \ — e) = 1, then there exists a 
nonreversible stationary distribution. 

Remark: Part (b) clearly follows from Corollary 2 if we replace the i.i.d. assumption with 
the assumption that {pi,i G Z 1 } be stationary and ergodic. The same is not true for part 
(a). To see this, consider the case of deterministic p^s with 



Then all stationary distributions are reversible if and only if a + (3 = 1. Indeed, if a + /3 = 1, 
then the process is symmetric, so all stationary distributions are exchangeable (and therefore 
reversible in this case) by Theorem 1.12 of Chapter VIII of Liggett (1985). On the other 
hand, the homogeneous product measures v p are stationary for all choices of a and j3 by 
Theorem 2.1(a) of the same chapter. However, since 



translate with probability ^ one obtains a stationary, ergodic sequence with nonreversible 
stationary distributions. So, in order to conclude that all stationary distributions are re- 
versible, it is not enough to assume that Pi > \ and pi < |, each a positive proportion of 
the time. 



(1.4) 




a if i is even, 
3 if i is odd. 
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2. Sufficient conditions for all stationary distributions to be reversible 

In this section, we prove Theorem 1 and Proposition 1. 

Proof of Theorem 1. The proof is based on the proofs of Theorems 1.2 and 1.3 of Liggett 
(1976), so we will only sketch the parts that are similar. We will use the basic coupling of 
two copies rit and r]' t of the process. In this coupling, particles move together as much as 
possible. At rate pi, for example, (77, rj') — > 

i) if v(i) = 1, r}(i + 1) = 0; = rf(i + 1) or rf(i) = 0, rf(i + 1) = 1, 

(»7, Vi,i+i) if rf{i) = 1, rf(i + 1) = 0; 77(1) = rj(i + 1) or 77(1) = 0, 7](i + 1) = 1, 

(»7M+i, < i+ i) if = ^(*) = 1; V(i + 1) = V(* + 1) = 0. 

Quantities related to the coupled process will be denoted by tildes. For example, the set of 
stationary distributions for the coupled process will be called X. 

For m < n, define the following two functions of a coupled configuration: 



fm, n M) = 1£ l \v(k)-rf(k)\, 



k=m 



9m,n(ViV') = # °f strict sign changes in the sequence {77' (m) — ry(m), ...,r}'(n) — rj(n)}. 

The fundamental property of these functions that makes them useful is that they cannot 
increase except as the result of transitions across the boundaries of [m, n] . Interior transitions 
can make them decrease unless rjt < rf' t or rjt > rj' t in the case of / m , n , and unless there is at 
most one sign change in the case of g m ,n- 

We will use the following notation. If v is a probability measure on {0, 1} Z x {0, 1} Z , 
then v{(r], rj') : rj(i) = S,rj'(i) = 5'} will be denoted by 




with analogous notation for probabilities of cylinder sets involving more than one site. 
If v e X, then J C,f m ^ n dv = 0. Writing this out gives 



2 J^(P< + 



Pm-l 



(2.1) 





+ similar terms coming from the right boundary of [m, n]. 
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The left side of (2.1) represents interior transitions that lead to the loss of (two) discrep- 
ancies, while the right side corresponds to the gain or loss of discrepancies in [m, n] due to 
transitions across the boundary. 

At this point, the argument differs somewhat according to whether one, both, or neither 
of the following hold: 



(2.2) 



inf(pi + q i+ i) > 0, m£(pi + q i+1 ) > 0. 

i<0 i>0 



If both of these infima are 0, then one can let m — > — oo, n — > +oo in (2.1) along appropriate 
subsequences to conclude that 



(2.3) 



(Pi + Qi+l) 



1 \ /0 1 
I 1 •(/ I 

i i + 1 / \i i+1 



= 0. 



It follows from this that v puts no mass on configurations with adjacent discrepancies of op- 
posite types, and then using the invariance of v again, that it puts no mass on configurations 
with discrepancies of opposite type at all. Therefore, it follows that 



(2.4) 



^{(77,7/) : r\ < rj' or 7/ < 77} = 1. 



We will now assume that (2.2) holds. (The argument in the case that one of the infima 
in (2.2) is zero and the other is positive is a combination of these two arguments, and will 
be omitted.) It now follows that 



00 




(I 


J ^ 




(0 




E 


V 









1 




i= — oo 






i + lj 






i + 1/ J 



< 00. 



Using the invariance of v again, it follows that 

E 



1 \ /0 
10 1 + v 1 

i i + k I \ i 




< 00 



for any k > 1. (See the proof of Lemma 4.4 of Liggett (1976) for details.) Next, using the 
fact that v e X implies j Cgm yn du = 0, one can show that v concentrates on configurations 
satisfying g m ,n(Vi v') ^ 1 f° r a ^ m < n. (See the proofs of Lemmas 4.7 and 4.8 of Liggett 
(1976).) 

The conclusion is the following: In all cases, v concentrates on configurations (77, rj') with 
the property that there is at most one strict sign change in the doubly infinite sequence 



{..., t/'(-2) - t/(-2), ri'(-l) - 7/(-l), 7/(0) - 77(0), r/(l) - 7/(1), t/(2) - 77(2), ...}, 



i.e., such that exactly one of the following is true: 
(a) 77 = 77', 
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(b) rf < T)' and r\ ^ rj', 

(c) rj >rj' and rj ^ rj', 

(d) there is a k so that rj(i) < rj'(i) for all i < k with infinitely many strict inequalities, 
and r](i) > rj' (i) for all i > k with infinitely many strict inequalities, 

(e) there is a k so that r/(i) > rj' (i) for all i < k with infinitely many strict inequalities, 
and 77(7) < rj'{i) for all i > k with infinitely many strict inequalities. 

(The fact that finitely many strict inequalities is excluded in cases (d) and (e) is a conse- 
quence of the fact that the system is in equilibrium; if there were finitely many, there would 
be some rate at which they would be annihilated, and since they cannot be created, this 
would contradict stationarity.) Since each of the above sets of configurations is closed for 
the evolution of the coupled process, if v e X e , then v concentrates on exactly one of these 
sets. 

Now take /j, fi' E T e such that (fi(fi) = (fi(fi') = 0. By Proposition 2.14 of Chapter VIII 
of Liggett (1985), there is a v e X e with marginals fx and fx' respectively. We will show 
next that v cannot concentrate on either of the last two sets described above. Suppose, for 
example, that it concentrates on the set described in (d). The k appearing there is random, 
so we will denote it by K(r),r)'). Let 

u m =p m - 1 fi{ri : r]{m - 1) = 1, 77(771) = 0} = q m n{ri : r]{m - 1) = 0, 77(771) = 1}, 
u' m =p m - 1 fu'{r] : 77(777 - 1) = 1, 77(777) = 0} = q m fi'{ri : r){m - 1) = 0, 77(777) = 1}. 

The second equality in each case comes from the fact that the fluxes are zero. Then 



Pm-1 



: 77(777 


"I) 


1 


1 


1 





777 — 1 


777 



\ / 1 



777 — 1 777/ \ 777 — 1 777 



1 \ / 1 \ / 1 
// I 1 1 iA - v\ lj. 

777 — 1 777 / y?77 — 1 777 J \ 771 — 1 777 

The third term on the right above is zero since v concentrates on the set described in (d). 
The second, fourth and sixth terms are bounded by 

(2.5) H(V,V') -K( V , rj') <m}. 

For the sixth term, for example, note that rj'(m — 1) = 1, 7/(777) = 0, 77(777 — 1) = 0, 77(777) = 1 
implies that K(rj, ?/) = m—1. The probability in (2.5) tends to zero as m — > —00. Therefore, 



Pm-1 





(u m - u' m ) -> 
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as m — > — oo. Similarly 





(u m - u' m ) -> 



as m — > — oo. Taking differences, we see that 



Pm-l 







tends to zero asm^ — oo. Note that this says that the sum of the first, second, fifth and 
sixth terms on the right side of (2.1) tends to zero as m — > — oo. The third, fourth, seventh 
and eighth terms are bounded by (2.5), so they tend to zero individually as m — > — oo. 
Applying a similar argument to the terms in (2.1) that come from the right boundary of 
[m, n], we see that the entire right side of (2.1) tends to zero as m — > — oo, n — > oo. It follows 
by passing to the limit in (2.1) that (2.3) holds, and then that (2.4) holds, contradicting the 
assumption that v concentrates on the set described in (d). A similar argument shows that 
(e) cannot hold either. Therefore, either \x < // or \x' < fi. 

To complete the proof of Theorem 1, take \x G X e and let \x' be one of the extremal 
reversible measures — v a or fi n according to whether (1.2) holds or not. For each such 
choice of //, we now know that either fi < fi' or fi' < fi. Suppose first that (1.2) holds. Then 
there is an a* G [0, 1] so that 



fi 



< v a if «o > ol* 
> v a if a < a*. 



Since the z/ a 's, together with the pointmasses on the configurations that are = and = 1, 
form a weakly continuous one parameter family, it follows that fi = v a * if a* G (0, 1), and 
is the pointmass on the = or = 1 configuration if a* = or 1. 

If (1.2) fails, then it follows that fi concentrates on U n A n U{0, 1}, where and 1 represent 
the identically and identically 1 configurations respectively, so that the statement of the 
theorem holds in this case as well. 



Proof of Proposition 1. We will prove the first statement; the proof of the second is similar. 
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Suppose that pi > \ — e for k < i < k + n, and take \i G X. Summing (1.3) gives 

k+n-l 

n(p(p) = [pMv ■ V(i) = hv(i + 1) = °) ~ Qi+ii^iv ■ v(i) = 0,77(7 + 1) = 1}] 

i=k 

> Q " e ) V{V ■ = 1, V(i + 1) = 0} - Q + : 77(1) = 0, r/(i + 1) = 1} 

k+n-l 

= -2e ^ ^{7/:7/(i) = l,7/(i+l) = 0} 

2+ e ) E [^{77 : 77(z) = 1} - yu{77 : 77(2+ 1) = 1}] 
' i=k 

> - 2en + Q + [^{77 : 7/(fc) = 1} - ^{77 : V (k + n) = 1}] . 

Since n is arbitrary, it follows that 4>(fS) > — 2e. To remove the extra factor of two, it suffices 
to note that 

k+n-l 

//{7/:7/(i) = l,7/(i + l) = 0}<^ 

for even n, since the events occuring in consecutive sumands above are disjoint. 

3. Sufficient conditions for the existence 
of nonreversible stationary distributions 

In this section, we prove Theorem 2 and Corollary 2. 

Proof of Theorem 2. Consider two choices p m -i, ■■■,p n and p' m _i, ■ ■■,p' n of jump probabilities 
satisfying p\ > pi for each i. Quantities related to the process corresponding to the p'^s will 
be identified with a prime. We need to show that 4>{^'m,n) — <fi(f J >m,n)- In order to do so, we 
construct a coupled process as follows: (X t ,r)t, rj' t , Y t ) has state space 

r n 

I (x, 77, r/, y) E Z 1 x {0, 1}^ x {0, 1}^ xZ 1 :x+^ [77(7) - 7/(1)] + y = 

k . 

and a; + ^ [//(i) — 7/'(i)] > for all m — 1 < k < n. > 

i=m 

Transitions inside [m, n] correspond to letting particles in the two configurations move to- 
gether as much as possible. Transitions across the boundaries (m— 1, m) and (n, n + 1) follow 
the same rules, with Xt and Yi keeping track of the number of discrepancies that leave or 
enter [m, n] at the left or right, respectively. To be more explicit, if the 1 appearing below 
are at sites i, % + 1 respectively (with m < i < n), then 

( f 1 \ 

• • • • • • at rate — p» 

v'\ _ ( 1 \ I V 1 J Pl 
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while if the 1 appearing below are at sites i, i + 1 respectively, for example, then 



/VW o i \ f (■■-! \") atrate ^-^i 

w v"o i"7 1 ^ i o \ 

••j o " ' J atrategj +1 . 

Note that in both cases, the discrepancy (^^j ^ s produced to the left of the discrepancy 

J^j . This means that the inequality x + J2i=m[v(^) ~ v'(^)] > is not violated by these 
transitions. 

At the left boundary, one has, for example, the following transitions: 

x ,o ...)) JHG •)) —a..-*- 



x, ( ••• J ] at rate p m _i, 



or 



X \VJ J = \ V '"/ / ~" V _1 ' V 1 " / J at rate P"»-l- 

While we have not listed all the possible transitions, hopefully we have listed enough so that 
the reader can easily construct the others. 

Now, start this process at (0, 77, 77, 0) , where 77 is any point in {0, l}[ m > n l. The limiting 
distribution as t — > 00 of r\t is fi m ,n while the limiting distribution of rj' t is fi' m n . For fixed 
m < k < n, let 

N t = (the number of times a particle in i]t has crossed from k to k + 1 by time t) 
— (the number of times a particle in rjt has crossed from k + 1 to k by time £), 

with iV/ being defined in an analogous way in terms of the process rj' t . Then one can easily 
check that 



Nl-N t = X t +J2lVt(i)-v't(i)}, 



so that N[ > N t a.s. On the other hand, 

j t EN t = p k P^[ Vt (k) = 1, rj t (k + 1) = 0] - q k+1 P^[ Vt (k) = 0, Vt (k + 1) = 1] 

so that 

d EN t 
<P{^m,n) = hni -J- EN t = lim — — . 

t^oo at t^oo t 

It follows from these observations that 4>{^' m n ) > 4>{,Hm,n) as required. 
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Proof of Corollary 2. Let \i m ,n be the stationary measure for the process on [m, n] corre- 
sponding to the given p^'s, and /j,' m n be the one for the process with p\ = \ + e. By Theorem 
2, <p(fim,n) > Wm,n)- By Theorem 2.9 of Liggett (1977), 



lim VL,n = "l/2, 

and therefore, 

lim ) = -• 

m^ — oo Z 

It follows that any weak limit [i of p m ,n as m — > — oo, n — > +oo is in X and satisfies > f . 
The measure is not reversible since all reversible measures have zero flux. 

4. The noninteracting case 

In this section, we consider a system of independent particles, each of which evolves as 
a continuous time Markov chain on Z 1 with unit exponential holding times and transition 
probabilities pi and qi from % to % + 1 and i — 1 respectively. It has been known since 
at least the publication of Doob's classic 1953 book that one way to construct stationary 
distributions for this system is to let {r](i), i £ Z 1 } be independent Poisson random variables 
with Erj(i) = <7j, where a = {o~i} is an invariant measure for the corresponding one-particle 
motion: 

(4.1) Gi = Oi-YPi-x + a i+1 q i+1 . 

By Theorem 4.12 of Liggett (1978), these provide all of the extremal stationary distributions 
if the one-particle chain is not positive recurrent. 

In such a stationary distribution, the flux of particles between i and i + 1 is 

(4.2) (p = piVi - q i+1 a i+1 , 

which is independent of i by (4.1). 
Solving (4.2) recursively leads to 

0~ n = CTm |H 1 h 



Qn - " - Qm+1 Qn \ Qn-1 Qn-1 ' ' ' Qm+1 

for m < n. It follows that if <j> > 0, there is a positive solution a to (4.2) if and only if 

1H 1 1 < oo, 

Po PoPi 

while if (j) < 0, there is a positive solution to (4.2) if and only if 

, , Po . PoP-i . 

1 H 1 1 < oo. 

9o 9o 9-i 



We therefore have the following result: 
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Theorem 4. Suppose the pi's are i.i.d. Then there is an extremal stationary distribution 
for the independent particle system with positive flux if Elog(po/qo) > 0, and one with 
negative flux if E \og(qo/po) > 0. 

Not surprisingly, these are exactly the conditions for a random walk in a random en- 
vironment to be transient to the right or left respectively - see Theorem 1.7 in Solomon 
(1975). 
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